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The generation of ultra-relativistic positron beams with short duration (t^+ < 30 fs), small 



divergence {0^+ ^ 3 mrad), and high density {n^+ 



10^ 



10 cm ) from a fully optical setup 



is reported. The detected positron beam propagates with a high-density electron beam and 7-rays 
of similar spectral shape and peak energy, thus closely resembling the structure of an astrophysical 
leptonic jet. It is envisaged that this experimental evidence, besides the intrinsic relevance to laser- 
driven particle acceleration, may open the pathway for the small-scale study of astrophysical leptonic 
jets in the laboratory. 



Creating and characterizing high-density beams of rel- 
ativistic positrons in the laboratory is of paramount im- 
portance in experimental physics, due to their direct ap- 
plication to a wide range of physical subjects, including 
nuclear physics, particle physics, and laboratory astro- 
physics. Arguably, the most practical way to generate 
them is to exploit the electromagnetic cascade initiated 
by the propagation of an ultra-relativistic electron beam 
through a high-Z solid. This process is exploited to gen- 
erate low-energy positrons in injector systems for con- 
ventional accelerators such as the Electron-Positron Col- 
lider (LEP) [1 . In this case, an ultra-relativistic elec- 
tron beam {E^- ^ 200 MeV) was pre- accelerated by a 
LIN AC and then directed onto a tungsten target. The 
resulting positron population, after due accumulation in 
a storage ring, was further accelerated by a conventional, 
large-scale {R ~ 27 km), synchrotron accelerator up to a 
peak energy of 209 GeV. The large cost and size of these 
machines have motivated the study of alternative particle 
accelerator schemes. A particularly compact and promis- 
ing system is represented by plasma devices which can 
support much higher accelerating fields (of the order of 
100s of GV/m, compared to MV/m in solid-state acceler- 
ators) and thus significantly shorten the overall size of the 
accelerator. Laser-driven generation of electron beams 
with energies per particle reaching [2-5 , and exceeding 
[6], 1 GeV have been experimentally demonstrated and 
the production of electron beams with energies approach- 
ing 100 GeV is envisaged for the next generation of high- 
power lasers (1-10 PW) [7J. Hybrid schemes have also 
been proposed and successfully tested in first proof-of- 
principle experiments [3 [9]. On the other hand, laser- 
driven low energy positrons {Eq+ ^1 — 5 MeV) have been 
first experimentally obtained by C. Gahn and coworkers 
[10 and recently generated during the interaction of a 
picosecond, kiloJoule class laser with thick gold targets 



pTHT4] . Despite the intrinsic interest of these results, 
the low energy and broad divergence reported {Eq+ < 20 
MeV and 6^+ > 350 mrad , respectively) still represent 
clear limitations for future use in hybrid machines. 

The possibility of generating high density and high en- 
ergy electron-positron beams is of central importance also 
for astrophysics, due to their similarity to jets of long 
gamma-ray bursts (GRBs), which are ejected as a conse- 
quence of the collapse of super-massive stars [15^ . These 
structures, despite extensive numerical studies [16], still 
present enigmatic features which are virtually impossi- 
ble to address by simply relying on direct observations. 
A possible solution might be represented by reproduc- 
ing small scale electron-positron jets (required bulk flow 
Lorentz factor of the order of 100 - 1000) in the labo- 
ratory. Although GRB jets may well have a weak large 
scale magnetic field [17 , the external shock is exclusively 
mediated by self-generated micro-scale magnetic fields. 
A purely electronic jet would present toroidal magnetic 
fields whose strength and structure would be compara- 
ble to the microscale fields that develop in response to 
the filamentation instability [18] and modify the shock 
physics. The presence of the highly mobile positrons 
would reduce the overall magnetisation of the jet (along- 
side with the amplitude of the electrostatic fields driven 
by a charge separation [19 ), simplifying the interpreta- 
tion of the experimental data and their comparison with 
the astrophysical scenario. 

Here we experimentally demonstrate the possibility 
of producing collimated and high-density ultrarelativis- 
tic positron beams in a fully laser-driven configura- 
tion. Their low divergence and short duration (com- 
parable to those of the incoming laser-driven electron 
beam) suggest the possibility of applying this generation 
scheme to current laser facilities towards the construc- 
tion of a fully-optical generator of high-quality, ultra- 
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Figure 1. a. Top-view of the experimental setup. Plastic and Teflon shielding was inserted to reduce the noise due to low 
energy divergent particles and x-rays. b. Typical positron signal as recorded by the Image Plate. The region labelled with 7 
noise is predominantly exposed by the 7-rays escaping the solid target, c. Typical signal of the electron beam as recorded on 
the LANEX screen, without a solid target and d. extracted spectrum. 



relativistic positron beams. Further more, the measured 
high positron Lorentz factors (7e+ c^ 200 — 400, com- 
pared to 7e+ < 10 in [10-12 ) of these beams are fi- 
nally comparable to those of astrophysical leptonic jets. 
This, in conjunction with the low divergence, the inferred 
electron/positron balance, and co-propagation with in- 
tense gamma-rays, finally open up a realistic possibility 
of studying the dynamics of such jets in the laboratory. 

The experiment (shown schematically in Fig. flla) was 
carried out using the HERCULES laser system at the 
Center for Ultrafast Optical Science (CUOS) in the Uni- 
versity of Michigan [20], which delivered a laser beam 
with a central wavelength Xl = 0.8 /im, energy El = 0.8 
J and duration tl = 30 fs. This laser beam was fo- 
cussed, using an //20 off-axis parabola, onto the edge of 
a 3 mm wide supersonic He gas-jet, doped with 2.5% of 
N2, with a backing pressure of 5.5 bar. Once fully ionised, 
this corresponds to an electron density of 9 x 10^^ cm~^. 
The focal spot size was measured to have a radius of 23 
/im which contained 50% of the laser energy (peak inten- 
sity of /l ~ 6 X 10^^ W/cm^). Laser power and gas-jet 
pressure were chosen in order to stay slightly above the 
threshold for ionisation injection [21 . This interaction 
delivered a reproducible electron beam with a divergence 
of approximately 2.5 mrad (see Figs, flld and file). Its 
spectrum was measured, prior to any shot with a high-Z 
solid target, by a magnetic spectrometer consisting of a 
0.8 T, 15 cm long pair of magnets and a LANEX screen. 
The arrangement of the spectrometer did not allow us to 
resolve electron energies below 80 MeV. Typical spectra, 
obtained using the calibration curves reported in [22], 
indicated the charge carried by electrons with energy ex- 
ceeding 80 MeV to be of the order of 50 pC (3 x 10^ 
electrons). Electron bunches obtained in similar condi- 
tions have been shown to have a length comparable to 
a plasma wavelength (Ape = ^TTc/ujpe ~ 10 jam) imply- 
ing a typical temporal duration comparable to that of 
the laser pulse [23 .The laser-accelerated electron beam 



interacted with mm-size high-Z solid targets of different 
materials (Cu, Sn, Ta, Pb) and thicknesses (from 1.4 to 
6.4 mm). The same magnetic spectrometer was used to 
separate the electrons from the positrons which were then 
recorded onto an Image Plate (IP). Due to the small dif- 
ference in positron and electron stopping power (below 
2% [24 ), the signal recorded was absolutely calibrated 
by using the calibration curves reported in [25 . Plas- 
tic shielding was inserted to reduce the noise on the IP 
induced by both the low-energy electrons and gamma- 
rays generated, at wide angles, during the laser-gas and 
electron-solid target interactions (see Fig. flla). 

In these experimental conditions, the positrons inside 
the high-Z target can be generated via either direct elec- 
troproduction (trident process), in which pair production 
is mediated by a virtual photon in the electron field [26] , 
or via a two-step "cascade" process where the electron 
first emits a real photon (bremsstrahlung) [27] which then 
produces an electron/positron pair via the Bethe-Heitler 
process [28 . Higher-order multi-step cascade processes 
may also significantly contribute to pair production, de- 
pending on the ratio between the thickness of the target 
and the radiation length of the material [29 . Keeping the 
parameters of the electron beam constant, the positron 
yield A/'g+ is expected to scale as: Nq+ cx {Z'^ndy ^ where 
n is the number of atoms per unit volume in the mate- 
rial, d is the thickness of the solid target, and j = 1 for 
the trident process and j = 2 for the two-step cascade 
process (we neglect here Coulomb corrections, which de- 
pend on Za, with a ^ 1/137 being the fine-structure 
constant). Neglecting the difference between the proton 
and the neutron mass, the mass density of the solid tar- 



get is p : 



AmpTi^ 



with A and nip being the atomic num- 



ber and the proton mass, respectively. If we maintain 
the areal mass density {a = pd) constant, we can then 
express the scaling as Np cx (Z'^/Ay . We have thus per- 
formed a series of shots for different materials (Cu, Sn, 
Ta, Pb) adjusting the target thickness so that the areal 
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Figure 2. Experimental (solid lines) and simulated (dashed 
lines) positron spectra for (a) 5.3 mm of Cu, (b) 6.4 mm of 
Sn, (c) 2.8 mm of Ta and (d) 4.1 mm of Pb. 



mass density was kept constant for each material (a « 4.7 
g/cm^, see first four rows in Table [l|). All the measured 
positron spectra, inferred from the signal recorded on the 
IP (which was absent if no solid target was inserted in 
the electron beam path, see Supplementary Material), 
showed a monotonically decreasing profile with approxi- 
mately 10^ positrons/MeV (solid lines in Fig. ^. In all 
cases, numerical simulations performed with the nuclear 
physics Monte-Carlo code FLUKA, which accounts for 
electromagnetic cascades during the passage of an elec- 
tron beam through a solid target [30] , are able to repro- 
duce the experimental data well (dashed lines in Fig. [2]). 
As theoretically predicted, the positron yield increases 
for materials with higher atomic number. This trend is 
quantitatively confirmed by integrating the experimental 
spectra in the range 90 < £;e+(MeV)< 120 (see Table 
[l] and Fig. |3|. Within this energy range, a maximum 
positron number of (2.30 ±0.28) x 10^ is obtained for the 
material with the highest Z (Pb). Keeping j as a free 
parameter, we obtain these data to be best fitted, as a 
function of Z^/A, if j = 2.1 ± 0.1 confirming the domi- 
nance of the cascade process with respect to the trident 
one (see Fig. [3Jb). The positron yield over the entire 
positron spectrum, as extracted from matching FLUKA 
simulations {Nt in Table [l]) , is seen to follow a similar 
trend. A further indication as to what process domi- 
nates is given by the dependence of the positron yield 
upon the target thickness {N^+ ex d for the trident and 
A^e+ ^ d^ foi" the two-step process) . A series of shots was 
thus taken varying the thickness of the solid target d for 
Ta and Pb (second four rows in Table IT] and Fig. Isla). 
As expected, the positron yield is seen to scale with d'^ 
in both cases. In order to support this statement the- 
oretically we compare, for each material, the radiation 
length Lrad with the range of target thicknesses d used 
in the experiment. The contribution of two-step process 
is expected to exceed the one of the trident process if 



<^/^rad ^ 10~^ ES]- For an order-of magnitude estimate 
of Lrad, we can assume here to be in the total-screening 
regime. For an electron with energy e emitting a photon 
with energy fkj, there is total screening if the parame- 
ter S = OLZ^I'^e{e — hhj)/ {hwmc^^ is much larger than 
unity, where the Thomas-Fermi model of the atom has 
been assumed [29]. In order to qualitatively explain the 
reported experimental results, we can roughly estimate 
e ^ huj ^ 100 MeV. Thus, it is 5* > 4 in the worst case of 
Cu, which is sufficiently large for the present estimate. In 
this regime, and by including Coulomb corrections, the 
radiation length is approximately given by [29] 



^rad — 



4a(Za)2nA2.Lo' 



(1) 



where n is the number of atoms per unit volume, 
Ac = h/mc = 3.9 X 10~^^ cm is the Compton wave- 
length, and Lo = log(183Z-i/3) - f{Za), with f{x) = 
Zl^i^^/^(^^ + x^). Eq. ^ provides the following 
values for the materials employed in the experiment: 
^rad(Cu) = 15 mm, Lrad(Sn) = 12 mm, Lrad(Ta) = 
4.1 mm and Lrad(Pb) = 5.6 mm. As seen by comparing 
these values with those in Table I, the material thick- 
nesses are always such that the inequality d/L^ad ^ 10~^ 
is fulfilled, in agreement with the experimental indication 
of the predominance of a two-step process for the electro- 
magnetic cascade. Moreover, in all the considered cases, 
except one where d = 4.2 mm for Ta, it is d < L^ad, 
which implies that the contribution of higher-order cas- 
cade processes can generally be neglected for an order- 
of-magnitude estimate. This is also corroborated by the 
observed angular divergence of the positron beams of the 
order of I/7, with 7 ^ 300 being the relativistic factor of 
the incoming electron beam (a larger number of cascade 
steps implies in general a broader angular distribution of 
the final particles). 

Due to the divergence of the positron beam, its max- 
imum density is located at the close vicinity of the rear 
side of the solid target. Here, the positron beam has a 



Mat. 



Cu 
Sn 
Ta 
Pb 



Ta 
Ta 
Pb 
Pb 



d(mm) 



5.3 
6.4 

2.8 
4.2 



1.4 
4.2 
2.2 

2.8 



9e+ (mrad) 
2.3 ±0.2 
2.7 ±0.3 
2.7 ±0.3 
3.5 ±0.4 



2.3 ±0.2 
2.7 ±0.3 
3.0 ±0.3 
3.3 ±0.3 



A^exp xlQ^ 



0.3 ±0.1 
0.6 ±0.1 
2.1 ±0.3 
2.3 ±0.3 



0.8 ±0.2 
3.8 ±0.3 
0.7 ±0.2 
1.1 ±0.3 



N • 



xlO^ 



0.3 
0.6 
2.1 
2.3 



0.8 
3.9 
0.7 
1.1 



Nt X 10^ 



31 
63 
190 
240 



78 
350 
60 
122 



Table I. The first four rows illustrate the results from targets 
with same areal density. The positron yield A^exp and A^gji^i 
refer to the energy window 90 < E^+ (MeV)< 120 as obtained 
from the experiment and FLUKA simulation, respectively. 
Nt refers instead to the total yield of positrons with E^+ > 
IMeV, as extracted from matching numerical simulations. 



transverse diameter of the order of 150 /im and, by as- 
suming that the positron beam wih retain the temporal 
duration of the initial electron beam {r^- < 30 fs, [23 ), 
a longitudinal length of the order of CTg- < 9 /im. In the 
case of maximum yield (4.2 mm Ta, see Table [l|) the den- 
sity of positrons with an energy between 90 and 120 MeV 
is of the order of 2.3 x 10^^ cm~^. FLUKA simulations 
indicate that this energy window contains approximately 
1% of the total positron yield. For 4.2 mm of Ta, this 
means that the total amount of positrons with energy 
Ee+ > 1 MeV win be of the order of 3.5 x 10^, indicating 
an overall positron density of about 2 x 10^^ cm~^. The 
overall positron beam intensity can thus be estimated 
to be of the order of 10^^ erg s~^cm~^. FLUKA sim- 
ulations show that such a positron beam co-propagates 
with an electron beam with an average density of about 
ne ~ 2 X 10^^ cm~^. The positron contribution on the 
leptonic beam will therefore be of the order of 10% with 
a null component of positive ions. 

We compare now our experimental results with the 
electron-positron astrophysical jets. Even though a de- 
bate is still open as to whether these jets are predom- 
inantly constituted by an electron-proton plasma or by 
electron-positron pairs, an element in favor to the latter 
is the power-law continuum spectra of the gamma-ray 
bursts associated with these structures (indication of the 
jets being optically thin) [31 . Despite the different gen- 
eration mechanism (pair production from gamma-gamma 
instead of gamma- nucleus collisions), their composition 
would be similar to the jets reported here, also thanks 
to their co-propagation with a high-density gamma-ray 
beam of similar size and duration (FLUKA simulations 
indicate a gamma-ray brilliance of the order of 10^^ — 10^^ 
ph/s/mm^/mrad^/0.1%BW). In the experimental results 
reported here, the excess of electrons in the beam implies 
a net current density of the order of Jg ^ —{'n~ —n'^)ec^ 
— 10''^"'^ A/m^ (assuming ng ~ 2 x 10^^ cm~^) inducing an 
azimuthal magnetic field of the order of l^^l ^ 30 T. 
However, FLUKA simulations indicate that, by varying 
the electron beam characteristics and target thickness, it 
is possible to significantly modify this percentage. For 
instance, the interaction of a GeV-like laser- accelerated 
electron beam with a thicker tantalum target {d ^ 2 
cm) is expected to generate a high-density, purely neu- 
tral electron-positron beam with an overall leptonic den- 
sity of the order of 10^^ — 10^^ cm~^. In this case, the 
toroidal magnetic fields would be virtually zero, allowing 
one to unveil the microphysics induced by the presence 
of small-scale magnetic fields generated by filamentation 
instability [11]. The high leptonic density would in fact 
allow for the study of the propagation of these jets in 
much more rarefied plasmas, in a scenario comparable to 
the propagation of these jets in the interstellar medium. 

The proposed mechanism for the generation of ultra- 
relativistic positron beams, applied to the near-term de- 
velopments in laser technology, might also be relevant 
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Figure 3. (a) Measured positron yield, in the energy range 
90 < E^+{Mey)< 120, for Ta (blue full circles) and Pb (red 
empty circles) for different target thicknesses. Lines give the 
best quadratic fits, (b) Measured positron yield, in the en- 
ergy range 90 < ^g+(MeV)< 120, for different materials but 
constant areal density as a function of Z^ j A. The dashed line 
represents the best quadratic fit. 



to the construction of all-optical electron-positron collid- 
ers. FLUKA simulations of the interaction of a pencil- 
like monoenergetic electron beam {E^- =100 GeV, over- 
all charge of 1 nC, see [7 ) with a 2 cm thick Ta target 
indicate the production of a positron beam with an ex- 
ponentially decreasing energy spectrum (maximum en- 
ergy of £^e+ = 80 GeV, with approximately 10^ positrons 
with energy between 70 and 80 GeV), a divergence of 
the order oi Q ^ 10 /irad and an overall charge compa- 
rable to that of the incoming electron beam. Particle- 
driven plasma wakefield accelerators can also be subse- 
quently employed to further increase the positron en- 
ergy in a metre-scale device, as recently demonstrated 
by Blumenfeld and collaborators [Sj. The extremely 
low- divergence achievable with our proposed generation 
mechanism would prove fundamental for efficient injec- 
tion of the positrons into such devices. The normalised 
emittance of such a positron beam can be expressed, for 
each monoenergetic component, as Sn ~ 7C^, being (" the 
beam source size. Due to the ultra-relativistic nature of 
the generation mechanism proposed, 7^ ^ 1 rad regard- 
less of the positron energy, thus reducing the normalised 
emittance to be 6^ ^ C ^ ^Ott mm mrad, in the conser- 
vative case of a 100 /im source size. This is comparable 
to the positron emittance measured after the injection 
stage of LEP {slep ^ SOtt mm mrad [33 ). It must also 
be noted that the positron beam would be inherently 
synchronised with the laser, allowing for the possibility 
of both electron-laser and positron-laser interactions. Di- 
rect comparison between these two experimental scenar- 
ios might allow for the testing of possible matter/anti- 
matter asymmetries in a highly non-linear regime. 
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